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I. INTRODUCTION 



Since the dawn of lattice field theory [j], the doubling problem of fermions has been a 
notorious obstacle to lattice simulations. Among several prescriptions for this problem, the 
Wilson fermion simply bypasses the no-go theorem 2| by introducing a species-splitting mass 
term into the naive lattice fermion. This Wilson term is regarded as one example of "flavored- 
maS8 te rmS » which split original 16 fermion species into plura. branch. flQ. It has been 
recently shown that the flavored-mass terms can also be constructed for staggered fermions 
[5-7] in Ref. js 10]. The original purpose of introducing these terms was establishment of 



the index theorem with staggered fermions |8j. A bonus here is that staggered fermions 
with the flavored-mass terms can be applied to lattice QCD simulations as Wilson fermion 
and an overlap kernel. One possible advantage of these novel formulations, called staggered- 
Wilson and staggered-overlap, is reduction of matrix sizes in the Dirac operators, which 



would lead to reduction of numerical costs in lattice simu 



ations. The numerical advantage 



in the staggered-overlap fermion have been shown in Now further study is required 

toward lattice QCD with these lattice fermions. 

The purpose of this work is to reveal properties of staggered- Wilson fermions in terms of 
the parity-phase structure (Aoki phase) 12H17I] . As is well-known, the existence of the Aoki 



phase and the second-order phase boundary in Wilson-type fermions enables us to perform 
lattice QCD simulations by taking a chiral limit since the critical behavior near the phase 
boundary reproduces massless pions and the PCAC relation. Besides, understanding the 



phase structure a 
and domain-wall 



so gives practical information for the application of the overlap 



20 



HQ 



2l| versions, both built on the Wilson-type kernel. Thus, in order to 



judge applicability of these new lattice fermions, it is essential to investigate the Aoki phase 
in the staggered- Wilson fermions. The phase structure for the staggered- Wilson fermion was 



22 



first studied by using the Gross-Neveu model in Ref. 
further investigation of this topic. 

In this paper, we investigate strong-coupling lattice QCD 
phase structure for two types of staggered- Wilson fermions 



231 ] and the present paper shows 







24|with emphasis on parity- 
101 ] . Firstly we discuss dis- 



crete symmetries of staggered- Wilson fermions, and show that physical parity and charge 
conjugation can be defined in both cases while hypercubic symmetry depends on types of 
staggered- Wilson fermions. Secondly, we perform hopping-parameter expansion and effec- 
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tive potential analysis for meson fields in the strong-coupling limit. For this purpose, we 
develop a method to derive the effective potential for lattice fermions with multiple-hopping 
terms. The gap equations show that pion condensate becomes non-zero in some range of a 
mass parameter, which indicates that parity-broken phase appears in this range. We also 
study meson masses around the second-order phase boundary, and find that massless pions 
and PCAC relations are reproduced. Lastly, we discuss parity-flavor symmetry breaking 
for 2-flavor cases. These results suggest that we can take a chiral limit by tuning a mass 
parameter in lattice QCD with staggered- Wilson fermions as with the Wilson fermion. 

This paper is organized as follows. In Sec. [TTJ, we review staggered flavored-mass terms 
and two types of staggered- Wilson fermions. In Sec. IIII1 we study discrete symmetries of 
staggered- Wilson fermions. In Sec. IIVI we study hopping parameter expansion in lattice 
QCD with these fermions. In Sec. [V] we investigate Aoki phase structure by effective 
potential analysis. In Sec. IVIt we discuss parity-flavor symmetry breaking in two-flavor 
cases. In Sec. IVIIj we devote ourselves to a summary and discussion. 

II. STAGGERED- WILSON FERMIONS 

Before looking into staggered- Wilson fermions, we review the Wilson fermion and its 
relatives. The Wilson term splits 16 species of naive fermions into 5 branches with 1, 4, 6, 
4 and 1 fermions. We call this kind of species-splitting terms "flavored-mass terms" . As 
shown in [jij], there are 4 types of flavored-mass terms for naive fermion which satisfy 75 
hermiticity. (75 in the naive fermion is flavored such as 75 £g> (T3 <8 T3 <S> r 3 <g> r 3 ) in the spin- 
flavor representation.) All these terms with proper mass shifts lead to a second derivative 
term as ~ a J dx^Dfa up to 0(a 2 ) errors. Thus we can regard them as cousins of Wilson 
fermion. 

There are also non-trivial flavored-mass terms for staggered fermions, which split 4 tastes 
into branches and satisfy 75 hermiticity. Since 75 is expressed in spin-taste representation as 
75 <8> 75 in this case, we only have two flavored-mass terms satisfying 75 hermiticity: 1 <S> 75 
and 1 (g> a^y. (For larger discrete symmetry one needs to take a proper sum for /i, v in the 
latter case.) These spin-flavor representations translate into four- and two-hopping terms in 



i 



the one-component staggered action up to 0(a) errors. The first type is given by 

M a = e Yl VimVBViCiC 2 C 3 C 4 = (1 <g> 75) + 0{a) , (1) 



syin 



with 



i?U = ( l) " 1 -- -'^.,, , (2) 
(%U = ( 1 1"'--"'- , (3) 
C^(V, + V^)/2, (4) 



The second type is given by 



M H = — = (rj 12 C 12 + VsaCu + r] 13 C 13 + r] 42 C 42 + VuC u + ^23^23) , 



1 

[1 <g> (0-12 + (T34 + o-i3 + o"42 + 0"i4 + 023)] + 0(a) , (6) 



with 



(7) 

(c^)^ = (-lr----'^ , (8) 
c p = (c,a+ac / ,)/2. (9) 



We refer to and Mh as the Adams- |8|] and Hoelbling-type 10[, respectively. The 
former splits the 4 tastes into two branches with positive mass and the other two with 
negative mass. These two branches correspond to +1 and —1 eigenvalues of 75 in the taste 
space. The latter splits them into one with positive mass, two with zero mass and the other 
one with negative mass. We note that Ma and Mh are also derived from the flavored mass 
terms for naive fermions through spin-diagonalization as shown in 3|. Now we introduce a 

n 

Wilson parameter r = r8 Xty and shift mass as in Wilson fermions [10[ . Then the Adams-type 
staggered- Wilson fermion action is given by 

S A = J2^[v^ + r(l + M A ) + M} xyXy , (10) 

xy 

D, = \(V,-Vi) . (11) 
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Here M stands for the usual taste-singlet mass (M = M5 x ^ y ). The Hoelbling-type fermion 
action is given by 

Sr = J2 Xx[v^ + r(2 + M H ) + M] xyXy . (12) 

It is obvious that these lattice fermions have possibility to be two- or one-flavor Wilson 
fermions. In lattice QCD simulation with these fermions, the mass parameter M will be 
tuned to take a chiral limit as Wilson fermions. For some negative value of the mass parame- 
ter: — 1 < M < for Adams-type and —2 < M < for Hoelbling-type, we obtain two-flavor 
and one-flavor overlap fermions respectively by using the overlap formula in principle. 

III. DISCRETE SYMMETRIES 

In this section, we discuss discrete symmetry of staggered- Wilson fermions. A potential 
problem with st agg ered- Wilson fermions in lattice QCD is discrete symmetry breaking. As 
discussed in 0, liol ] . the discrete symmetry possessed by the original staggered fermion is 
broken to its subgroup both in the Adams-type and Hoelbling-type actions. We below list 
all the staggered discrete symmetries (shift, axis reversal, rotation and charge conjugation), 
and look into their status in staggered- Wilson fermions. Shift transformation is given by 

3u : Xx -> (n(x)Xx+ii, Xx -> (p{x)Xx+fi, Uu,x -> U UtX+[l , (13) 

with Ci(ar) = (_i)*a+*s-He4 j = (-1)^+^ £ 3 ( x ) = (-1)** and ( A (x) = 1. It is obvious 

that this transformation flips the sign of both flavored-mass terms. The Adams type is 
invariant under the two-shift subgroup asx— > x + l±£i while the Hoelbling type is invariant 
under four-shift subgroup asx— )-x + i±2±3±4. Note that these subgroups include the 
doubled shift x — > x + 2fi as their subgroup. The axis reversal transformation is given by, 

Z» ■ Xx (-l) x "xix, Xx (-l) x "Xi x , Uu,x -»• U„, Ix , (14) 

with I = I 11 is the axis reversal x^ — > —x^, x p — >• x p , p ^ fi. It again flips the sign of the 
flavored-mass terms. The staggered rotational transformation is given by 

K pcT : Xx ->■ S r (R~ 1 x)xr-i x , Xx -> S^Rr^XBrix, U U)X ->■ U U)Rx , (15) 
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where R pa is the rotation x p — > x a , x a — > — x p , x T — > x T , r ^ p,a and Sr(x) = |[1 ± 
Vp( x )v<t( x ) T (p( x )(o-( x ) + ? 7p(^) r ?o-( a; )Cp( x )Co-( :;c )] with p ^ cr. It is notable that the Adams- 
type fermion keeps this staggered rotation symmetry while the Hoelbling type loses it. The 
staggered charge conjugation transformation is given by 

C : Xx -> e^xl, Xx -> -e^xl, -> ^ • (16) 

The Adams-type fermion keeps this symmetry while the Hoelbling type loses it. 

We next elucidate residual subgroups possessed by staggered- Wilson fermions, and discuss 
how to define physical discrete symmetries as Parity, Charge conjugation and Hypercubic 
symmetry. For this purpose we separate spin and flavor rotations in the above transforma- 



tions. Here we utilize the momentum-space representation in 25|, |26[. In this representation 
we can define two set of Clifford generators T p and S M , which operate on spinor and flavor 
spaces in the momentum-space field 4>(p), respectively. (Details are shown in Appendix. |A~1 ) 
Then the shift transformation translates into 

S M : <t>(p) -> exp(ip p )E p (f)(p) . (17) 
The axis reversal translates into 

x m :0(p) -> r M r 5 s M s 5 0(/ P ) . (is) 

The rotational transformation translates into 

TZ pa : <j>{p) -)• exp(^r p r (T )exp(^S p S (T )0( J R-V) • (19) 

By using this representation we figure out residual discrete symmetries of staggered- Wilson 
fermions as follows. We first consider parity. Both staggered- Wilson fermions are invariant 
under 

x s Si ~ exp(ip 4 )rir2r 3 r5</>(-p,p4) ~ exp(ip 4 )r 4 0(-p,p 4 ) , (20) 

with X s = X1X2X3. This is essentially the continuum parity transformation 



In the contin- 



uum limit the phase factor disappears and it results in the continuum parity transformations 
P : ij}(p) —> 7 4 ^(— P,p 4 ) for the Dirac fermion. We thus conclude both staggered- Wilson 
fermions possess symmetry leading to physical parity symmetry P. We note the simple 
combination of //-shift and /i-axis reversal (shifted-axis reversal :S p X p ) is also a symmetry of 
both fermions. 



7 



We next consider physical charge conjugation. In the case of Adams fermion the stag- 
gered charge conjugation symmetry C in Eq. ffl6|) remains intact. Thus, physical charge 
conjugation for the two-flavor branch can be formed in a similar way to usual staggered 



fermions as shown in |26[] (C ~ C1S21S4X2X4). On the other hand, the Hoelbling type breaks 
C. In this case, however, we can define another charge conjugation by combining C with 
rotation transformation as 

c T - n 21 nj 3 c . (21) 

The Hoelbling action is invariant under this transformation. By using this symmetry we can 
define physical charge conjugation C for one-flavor branch as in the Adams type. We thus 
conclude that both staggered- Wilson fermions have proper charge conjugation symmetry. 

TABLE I: Symmetries of staggered- Wilson fermions. 





N f 


5&X-subgroup 


7£-subgroup 


P C 


sw 4 


Staggered 


4 


Sfj,, If_t 




OO 


O 


Adams 


2 






OO 


O 


Boelbling 


1 






OO 


X 



We lastly consider hypercubic symmetry. In staggered 



'ermions, the rotation Eq. f JT5|) 
and axis reversal Eq. (fl4j) form hypercubic symmetry 27J, which enhances to Euclidian 
Lorentz symmetry in the continuum limit. In the case of the Adams-type fermion, the ac- 
tion is invariant under the rotation Eq. (fl3]) and the shifted- axis reversal S^X^. These two 
symmetries can form proper hypercubic symmetry SW4 in this case. Thus we conclude that 
Adams fermion recovers Lorentz symmetry in the continuum. In the Hoelbling-type formu- 
ation, the action breaks the rotation symmetry into its subgroup called doubled rotation 
lo| . which is given by 

n pa n, u ~ ex P [^(r p r CT + r^)] ex P [^(s p H CT + z^mR^R^p) , (22) 

where (fi,u,a,p) is any permutation of (1,2,3,4). It is also invariant under sequential 
transformations of (//, v rotation), (z/, fi rotation), (yu shift) and {y shift) as 

SvS^ll^K^ ~ exp(zp p + ip^T^Ty <j>{p) , (23) 

with = —Pp, u ,pT = Pt, t 7^ fi, v. The loss of rotation symmetry indicates that we 
cannot define hypercubic symmetry in the Hoelbling fermion. It implies that it could not 



s 



lead to a correct continuum theory, and we would need to tune parameters to restore Lorentz 
symmetry. Indeed the recent study on symmetries of staggered- Wilson fermions by Sharpe 
reports that recovery of Lorentz symmetry requires fine-tuning of coefficients in the 



gluonic sector in lattice QCD with Hoelbling fermion. 

To summarize, Adams fermion possesses physical parity, charge conjugation and hypercu- 
bic symmetries while Hoelbling fermion loses hypercubic as shown in Table. [H It seems that 
Hoelbling fermion cannot be straightforwardly applied to lattice QCD while Adams type 
can. We note that both staggered- Wilson fermions have proper parity symmetry, and we 
can discuss spontaneous parity symmetry breaking. Moreover, we may find some symptom 
due to Lorentz symmetry breaking in Hoelbling fermion in the parity-phase structure. This 
is another motivation to study parity-phase structure in strong-coupling lattice QCD. 

In the end of this section, we comment on special symmetries in staggered- Wilson fermions 
without mass shift. Hoelbling fermion without the mass shift (r]^D^ + M H ) possesses special 
charge conjugation symmetry {C' T : x ~^ Xi X ~ > x)- This topic is beyond the scope of this 
work, but we note that it can do good to two flavors in the central branch. Use of the central 



branch is intensively discussed 



in [n 



291. 



IV. HOPPING PARAMETER EXPANSION 

In this section we investigate parity-phase structure in lattice QCD with staggered- Wilson 
fermions in the framework of hopping parameter expansion (HPE) in the strong-coupling 



regime 12] . In the hopping parameter expansion, we treat a mass term as a leading action 
while we perturbatively treat hopping terms. We thus perform expansion by a hopping 
parameter which essentially corresponds to inverse of a mass parameter. By using self- 
consistent equations we derive one- and two-point functions, and calculate meson conden- 
sates and meson mass for two types of staggered- Wilson fermions. We for simplicity drop 
the flavor indices until we discuss the two-flavor case in details in Sec. IVII However it is easy 
to recover the flavor indices for the field Xfi the mass parameter Mf and the condensate £/ 
(/=1,2,...). 
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— hopping (mass term) 



(x,a) (y,b) 
1 — hopping (kinetic term) 



(x,a) (x + p,, b) 



(x,b) (x + fl,a) 
2 — hopping (flavored mass term) 



(X a x X b y)° = "M" 

Kr^ x (U^ x ) ab 

;(2) 



FIG. 1: Feynman rules in hopping parameter expansion (HPE) with the Hoelbling-type staggered- 



Wilson fermion. a and b stand for the color indices. VVl , ° is given in Table ITT1 



TABLE II: Concrete forms of W^, in Fig. ffl 



- - C/t .tft _ 

^J,,X — fl U,X — fl — U 

+ - u ux u ] - 

t 1 ^ v,x+p,—v 
+ U\,x-(Pv,x-(l 

A. Hoelbling type 

We begin with the Hoelbling-type fermion, which contains two-hopping terms in the 
action. One reason that we start with Hoelbling type regardless of the lower rotation sym- 
metry is that 2-hopping calculation can be a good exercise for 4-hopping case in Adams 
type. To perform the HPE for the Hoelbling-type fermion, we rewrite the action Eq. ( !T2|) 
by redefining x ~~ V%Kx with the hopping parameter K = 1/[2(M + 2r)], 

S = Yl * xXx + 2K XxiVuDjxyXy + 2Kr Xx(M H ) xy Xy , (24) 

x x,y x,y 

where Mh is given by Eq. (|6]). The plaquette action is 1/g 2 term and we can omit it in 
the strong-coupling limit. In this section, we derive one- and two-point functions by using a 
0(K 3 ) self-consistent equation: Solving this equation leads to truncation of diagrams as the 



10 




FIG. 2: Feynman diagram for mesonic one-point functions in the 0(K 3 ) self-consistent equation 
of HPE with the Hoelbling fermion. Black circles stand for the leading one-point function {xxXx)o 
while white circles stand for (XxXx) which include next-leading and higher hopping terms. By 
summing up higher contributions, we obtain the second equality. 

ladder approximation having all diagrams to 0(K 3 ) are taken into account. More precisely, 
this approximation does not take account of all diagrams, but it successfully includes certain 
kinds of diagrams to all orders of K thanks to a self-consistent approach. We thus expect that 
it works to figure out existence of Aoki phase. We note that this approximation especially 
works well for a small hopping parameter K <C 1. In Fig. [IJ we depict Feynman rules 
in the HPE for this fermion. The fundamental Feynman rules contain contributions from 
0-hopping (mass term), 1-hopping (kinetic term) and 2-hopping (flavored-mass term) terms. 

First, by using these Feynman rules, we derive meson condensates from an one-point 
function of the meson operator Ai x = XxXx m the mean-field approximation. The one-point 
function is defined as 

{XxXx)- <W^_ J v[x ^ u]eS ■ W 

Note that we use the partition function Z = f V[x, x, U]e s , not Z = f T>[x, X, U}e~ s , 



24|. The 



following the convention for the partition function in the strong-coupling analysis 
leading term in the hopping parameter expansion is given by 

\XxXx/0 - fri ,, - rn ~ ~° , \. zl °) 



fV[ X ,X,U]e 
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where So = J2 X XxXx- By using the Feynman rules, we can evaluate the diagrams in Fig. |2j 

(x a x xl) = -* ab Zx 
= (xlx% 

+ ^(-l)(^) 2 (( X a x)x)o^((xx), + A)o^((xX 6 ),)o 
±*» 

+ (- 1 ) ( 2 ^ r ^^7f) ((^)-)o^it((xx), + A + p)oW^(( X x b ) :[ .) 

+ ^ (-l)(^J 2 (-l)(AX,) 2 ((x a x),)o^((xx)x +A )o^((x^ 
x ((xx)x-+ ! ))o^i :E ((xX 6 )x-)o 

+ ± ( 2ifri ^^ ! ) 2( - 1) ( 2A '""-i^i) 2 

x ((x a x).)oW^ x .((x^ + ^ 

+ Yl (-^(^w^V) (-i)(%,) 2 ((fx),}o>vg, 

±/i,±i/,±p ^ V / 

x ((n)i+A#)o^p,i+(i+s(te)i+M+f+p)o^ + fi + ii ((xx)x+A +i >)o>V^((xX b )x)o 

+ 0(K 5 ) , (27) 

where (xx)x stands for XxXx and = W+^+j,^ in Table. [TTJ Note that we consider only 

connected diagrams in Fig. [2j By summing higher hopping terms, the one-point function is 
obtained as shown in Fig. [2, which is given by 

-E X = -(M X ) = -(M x ) + 2K 2 J2ZxZ x+ t-2-^ {Kr) 2 J2 E * E *+(>+fi ■ ( 28 ) 

The equation contains terms to OiK 2 \ and 0(K 3 ) diagrams are found to vanish due to 
cancellation between the diagrams. Here we solve it in a self-consistent way for condensate 
S within mean-field approximation. We here assume T, x = a x + ie x ir x as the condensate, a 
and 7r correspond to chiral and pion condensates, respectively. We substitute this form of 
S x in Eq. ( I2"gj) and obtain a self-consistent equation 

- (a + ie x ir) = -l + 2K 2 -4 (a 2 + tt 2 ) - 2 • ^{Kr) 2 • 4 • 3 (a + ie x n) 2 , (29) 
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' x * * * x " * 



fj, V 



fi + v 



FIG. 3: Feynman diagram for mesonic two-point functions for 0(K 3 ) self-consistent equation with 
the Hoelbling fermion. 

which yields —a = — 1 + lQK 2 n 2 and —in = —8K 2 ■ 2ian. For simplicity, we have set 
r = 2\/2 to make the equation Eq. (1291) simpler. Of course we can also discuss for other 
values of r in general. 

Now we have two solutions depending on n = or n / 0: For ir = 0, we have a trivial 
solution cr = 1. For it ^ 0, we have a solution as 



1 



a 



1QK 2 ' 



7T 



± 



/ 1 


V 1 16K 2 ) 


/ 1QK 2 





(30) 



Nonzero pion condensate implies spontaneous parity breaking for the range | K | > 1/4. The 
sign of the pion condensate in Eq. ( 1301) reflects the parity symmetry of the theory. Thus 
the parity-broken phase, if it exists, appears in a parameter range -4V2-2 < M < -Ay/2+2 
in the strong-coupling limit. We note that the critical hopping parameter \K C \ = 1/4 is small, 
and we speculate that the 0(K 3 ) self-consistent equation is valid around the value. 

Next, we discuss the meson mass from a two-point function of the meson operator 
S{0,x) = (AioAi x ). From Fig. [21 we derive the following 0(K 3 ) equation for a two-point 
function. 

S(0,x)=(x a XoX b x X b x ) 
= - S 0x N c 

-K 2 (x a oXoXo(v,,o) 2 



XoXxX x) 



a,/3=± [i^u 



O^a^ufi) Xafil3C>Xajll3u{YVa[lpv,0 



Xq XxXx) 



(31) 



13 



where Walavo ls defined in Table. [TTJ Note that we consider only connected diagrams in 
Fig. [3j By integrating out the link variables in the strong-coupling limit, it is simplified as 

5(0, x) = (x a XoX b x X b x ) = S 0x N c + K 2 X)0$®&£) 

+ ( 2 #™^=) Efe^>- ( 32 ) 
Then the self-consistent equation for S is given in the momentum space as 



S{p) 



-Nr. 



) 



+ (W r— ^ V (e-^^+f^ + 







5(p) 



(33) 



5(p) = iV c 



(34) 



We finally obtain the meson propagator as 

/ 1 \ 2 1 _1 

-2if 2 cosp M + 4 I 2A> - J ^cosp^cosfv-l 
M V 2 V3/ ^ 

Here the pole of S (p) should give meson mass. Since x is an one-component fermion, it may 
seem to be difficult to find the pion excitation from Eq. ( I34j) . However, as we discussed, 
75 in the staggered fermion is given by e x = (— \^ Xl+ - +X4 and the pion operator is given 
by tt x = Xxi^xXx- We therefore identify momentum of pion by measuring it from a shifted 
origin p = (it, it, it, it). Here we set p = (im n a + n, it, it, it) for 1/S(p) = in Eq. (|34|) . Then 
we derive the pion mass m n as 

1 - 16K 2 



cosh(m vr a) 



QK 2 



(35) 



where we again set r = 2a/2 for simplicity. In this result, the pion mass becomes zero 
at | A" | = 1/4, and tachyonic in the range | K |> 1/4. It implies that there occurs a 
phase transition between parity-symmetric and parity-broken phases at \K\ = 1/4, which is 
consistent with the result from the one-point function in Eq. (1301) . We note that the massless 
pion at the phase boundary is consistent with the scenario of second-order transition. We 
can also derive the sigma meson mass by substituting p = (im n a, 0, 0, 0) for 1/S(p) — in 
Eq. (EH) as 

1 



cosh(m CT a) 



1 + 



2K 2 



(36) 
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— hopping (mass term) 



(x,a) (y,b) 
1 — hopping (kinetic term) 



(x, a) (x + fi,b) 



(x,b) (x + p,a) 
4 — hopping (flavored mass term) 



lx,a) n 



a (x + ap + (30 + ~/p + 5a, b) 



(x,b) /Li 



p a (x + ap + (30 + 7 p + 5a, a) 



-Kv,Aul x ) ab 



FIG. 4: Feynman rules for the HPE with the Adams fermion. a and b stand for the color indices. 
We show the concrete forms of Sa x in Table IIIII of Appendix [Bj 

B. Adams type 

We investigate the parity-phase structure for the Adams-type staggered-fermion by using 
0(K 3 ) self-consistent equations in hopping parameter expansion. The approach is basically 
parallel to that of Hoelbling type. We just need to consider Feynman diagrams for this case. 
The action Eq. (fTUj) is rewritten by redefining x ~ ► V^Kx with K = 1/[2(M + r)] as, 

S = * xXx + 2K H Xx(Vn D n)*yXy + 2A> Xx(M A ) xy Xy , (37) 

x x,y x,y 

where Ma is given in Eq. (JT]). In Fig. HI the Feynman rules in the HPE for this fermion are 
depicted. First, we derive meson condensates from the one-point function hA. x — XxXx- 

The 

equation for the one-point function is obtained as shown in Fig. |SJ 
-E X = -(M X ) 

= -(M x ) + 2K 2 - 2 ■ 1 - 23 (Kr) 2 ]T E x E x+fl+i>+ ^ a . (38) 

We substitute H x = a x + ie x n x for E x in Eq. (138]) and obtain the self-consistent equation 
- (a + ie x ii) = -1 + 2K 2 ■ 4 (a 2 + tt 2 ) - 2 • — i-^ (Kr) 2 • 4! (a + ze^vr) 2 . (39) 

From this, we obtain —a = — l + lGA^ 2 ^ 2 and —in = —8K 2 -2iair. Here we have set r = 16a/3 
to make the equation simple. We again have two solutions: For n = 0, we have a trivial 
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FIG. 5: Feynman diagram for mesonic one-point functions in the O(K^) self-consistent equations 
of HPE with the Adams fermion. There is a 4-hopping fundamental diagram, which is peculiar to 
this fermion. 



' x " * * x * * 



u fi + v + p + a x 

FIG. 6: Feynman diagram for mesonic two-point functions for O(K^) self-consistent equations of 
HPE with the Adams fermion. 



solution a — 1. For 7r ^ 0, we have a non-trivial solution as 



a 



16K 2 



71 = ±n 



(40) 



It indicates that parity-broken phase appears in the range of the hopping parameter as 
| K |> 1/4 or equivalently -16^-2 < M < -16^+2. 

Next, we discuss the meson mass from a two-point function of the meson operator 
S(0,x) = (M.qM. x ). From Fig. [6] we derive the following equation for two-point functions, 
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S(0,x)=(x a X a oXlxl) 
= - 5 0x N c 

-K 2 (x a X a X c (v,,o) 2 
( 



u%xU e M,o) ef + K- fi rx%x e ^u;'_ 



Xo XxXxi 



;(4) 



cd 



X W 



(4)t 



+ ( 2-ft>e7/5 j (XoXoXo ^ ^ y^a^v^pSafl J Xap,/3v'yp5aXap,pi>'yp5a 

xfoM) , (41) 

where W^g fa ^ is defined in Table. llHl of Appendix [B] By integrating out the link variables 
in the strong-coupling limit, it is simplified as 

5(0, X ) = (Toxztxi) = - So*n c + k 2 J2(xinxixi) 



±n 

2 



(^^ r A \ . 9 4 ) (XjX+i>+p+aXfl+u+p+&XxXx) ■ (42) 

/ + ./ -i-n 4-rr 



Then the self-consistent equation for S is given in the momentum space as 



S{p) = -N c 



+ 1 2Kr- 



-K 2 ( e ~* PM + eiPM ) 



V 4! • 2 4 
We finally obtain the meson propagator as 



E E 

a,/3,7,(5=± n^iyjtp^cr 



(43) 



S(p) = iV c 



-2K 2 cosp^ + 16 yiKr—^- 



4! .24 



cos p M cos p v COS cos p a — 1 



(44) 



Here we set p = (im^a + 7T, 7r,7r,7r) for l/5(p) = in Eq. (|4*4"j) . which gives the pion mass 
m n as 



cosh (m n a) 



1 - 16K 2 
10K 2 



(45) 



where we again set r = 16 a/3 for simplicity. Here the pion mass becomes zero at | K \= 1/4 
and becomes tachyonic in the range | K |> 1/4. It suggests that there occurs a second- 
order phase transition between parity-symmetric and broken phases at \K\ = 1/4, which 
is consistent with Eq. f|4*0|) . We can also derive the sigma meson mass by substituting 
p = [im n a, 0, 0, 0) for 1/S(p) = in Eq. flH| as 

1 



cosh(m cr a) = 1 + 



6K 2 



(46) 
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V. EFFECTIVE POTENTIAL ANALYSIS 



In the previous section, we have investigated the parity-phase structure in hopping pa- 
rameter expansion. We found a strong sign of parity-broken phase for | K \ > 1/4. In order 
to judge whether the parity-broken phase is realized as a vacuum, the analysis of the gap 
solution in the hopping parameter expansion is not enough, and we need to investigate the 
effective potential for meson fields. 

In this section, we consider the effective potential of meson fields for SU(N) lattice gauge 
theory with staggered- Wilson fermions. In the strong-coupling limit and the large N limit, 
effective action can be exactly derived by integrating the link variables 12|, I24J- Then, by 
solving a saddle-point equation, we can investigate a vacuum and find meson condensations. 
In this section we again begin with the Hoelbling case as exercise, and go on to Adams 
fermion with better discrete symmetry. 



A. Hoelbling type 



In the strong-coupling limit we can drop plaquette action. Then the partition function 
for meson fields M. x = {XxXx)/N with the source J x is given by 



Z{J)= V[XiX,U}ex V 



N^J X M X + S F 



(47) 



where Sf stands for the fermion action. Here we have defined Ai x with a 1/N factor to 
ensure it to have a certain large N limit. In this case, Sf is the Hoelbling-type staggered- 
Wilson action Eq. ( f!2|) . N stands for the number of color. In the large ./V limit, we can 
perform the link integral. We here consider the effective action up to (9(A4 3 ) for meson field 
M.. This order corresponds to the 0(K 3 ) self-consistent equation in the hopping parameter 
expansion. 

We develop a method to perform the link-variable integral with multi-hopping fermion 
action terms. In our method, we perform the link integral by introducing two kinds of 
link-variable measures. Now we formally rewrite the partition function as, 



Z{J)= /P[ X ,x]exp 



Y^N^Jx + m) M Q 



cxp 



^AW(A) 



(48) 
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{Uv,x+p,-> U v x+ p) 



X fX + [l 



t 



FIG. 7: Link variables corresponding to the two kinds of measures in the partition function Eq. 
19 D in a 2 dimensional case. 



where we define M = M + 2r and the last term is expressed as, 
exp J]AW(A) =JI Z *' 

a: J a; 

= J ^JJ 2? [U^, U^ x+i> ]j exp [- (Tr(^t) - Tr(V^E))] . (49) 

A, E and £^ (V and V^) are composites of the fermion field x (link variables U), which 
we will explicitly show later. W(A) is a function of A, which will be an essential part of 
effective potential of meson fields. 

Now we explain how the integral in Eq. fj4*9]) can be performed by using two types of the 
link measure. Let us consider two-dimensional cases in Fig. [7]fbr simplicity. In this case, U^ x 
and Upx+o (fi v) form link variables in a square block. We can classify diagrams to 0(Ai 3 ) 
into three types: (1) 1-link (//) + 1-link (— jj) hoppings, (2) 2-link (//, v) + 2-link (— u, —fi) 
hoppings, (3) 2-link (/i, v) + 1-link {—v) + 1-link (—fi) hoppings. The 1-link hopping comes 
from the usual staggered kinetic term while the 2-link hopping from the flavored-mass term. 
(1) and (2) are 0(M. 2 ) while (3) is 0(A4 3 ). Since one square block contains all the three 
diagrams, we can derive the effective potential by integrating link variables per each block. 
We note that 0(Ai 3 ) diagrams cancel between one another, which is consistent with the 
HPE calculations. We can also avoid double-counting by adjusting factors for one- link and 
two-link terms as we will show in Eq. (I64j) . 

In this method, we need to define sets of link variables and fermion bilinears as V and 
E in Eq. f )49|) : V and E are matrices including components corresponding to 1- and 2-link 
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terms. We call a space spanned by these matrices "hopping space". Here we define a, b and 
a, /3 as color and hopping space indices respectively. We also denote Tr as trace for color 
and the hopping space. Explicit forms of V and E are given by 

V^ = diag{V^,V 2 ab ,Vf) , (50) 

with 

Vi = diag (U^ x ) 

= diag(U lx ,U 2x ,--- ,U 4>X ) , (51) 

V v ' 

1 

V 2 = diag {U^ x U VtX+[ i) 

= diag (U ljX U 2 +h U liX U 3 +h ■■■ , U 4iX U 3x+i ) , (52) 

V «, ' 

12 

V 3 = diag [U^+pUl^ 

- diag (r,,.. 2 r 2 ,..,. u 1>x+ji ul x+v r . ;; rj . ,) , (53) 

V v ' 

12 

E% = diag «, Ef, Ef) , (54) 



and 



A = diag (A, M ) 

= diag(A,i,A,2,--- ,A, 4 ) , (55) 



1 



A = diag (A,/m>) 

= diag (A,i2, A,i3, • • • , A,43), = 2, 3) , (56) 
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where we define the operator D as the fermion bilinears, 



\ ao 



ah 1 1 

1 „ h , n snh J- 

A,/^) ~ ^^y^VfJ.u,xX x X x +ll+u i (A (58) 
= ^^y^V^,x+uX x +uX x +p, i (A,/w) = 23^/g ? ^A tl/ I : ' ; +! > ^2;+/iXx+P ' (59) 
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Here V\ and E\ are 4x4 diagonal matrices while Vj, and Ei (i = 2, 3) are 12 x 12 diagonal 
matrices. Now we have prepared to obtain W(A). By using the relation U'U = 1, we obtain 
the Schwinger-Dyson equation, 

d 2 Z T 



be 



9E%d{E\) Pi 



(60) 



W should be a function of a gauge-invariant quantities as follows. 

A- = ± (E<eC ■ (61) 
We can solve the Schwinger-Dyson equation analytically and derive W as a function of A, 

^l + (l-4A) 1/2 ~ 



W(A) = Tr 



1_4A) 1/2 - 1 - In 



(62) 



We here perform trace for colors and hopping spaces. 



5X A ) = -E 



(l-4A x .) 1/2 -l-ln 



1 + (1 - 4A, 



4/2 



(63) 



Finally we obtain a concrete form of A as 
1 



A 



M x M x+il + \Y1 M *+t M 



x+[i+0 



L H 



—7= J] (MX + M x+i >M xH 

(64) 



The first and second terms correspond to contribution from Z>i )At , Z?| , and the third and 
forth terms correspond to contribution from D i fJi , D\ (i = 2, 3). Now we again set r = 2^2 
to compare the result to that of the hopping parameter expansion in Sec. HVl We need to 
change the fermion measure to the meson-field measure as 



Jv\x,x] = J VM exp 



(65) 



Then the effective partition function for the meson field is given by 



Z(J) = / DM exp 



Nfcj x M x + Srf{M)\ 



S cS (M) = (MM X - \nM x ) + W(A) 



(66) 
(67) 
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where we denote M as the shifted mass parameter M = M + 2r. The partition function 
with J = in the large N limit is reduced to the integrant for the saddle-point values of the 
meson fields, 



(68) 



Z(J = 0) = J VMexp [NS cS (M)] 

~ exp [NS cS (M)] , (N -> oo) . 



Now we consider pion condensate. For now, we consider only scalar a and pseudo-scalar 7r 
fields as 



M x = a + ie x ir , 



(69) 
(70) 



By substituting this form of the meson field into the Eq. ( IBTl) . we derive the effective action 
for the £ and 6, 



S cS (M) = M^£cos# ~ 5Z lnS 



2- 2 /A 1/2 



(l - 4-2E 2 sin 2 0) 



In 



2 /A 1/2 



1+ (l-4-2E 2 sin 2 0) 



(71) 



We ignore the irrelevant constant. From the translational invariance, we factorize the 4- 
dimensional volume V4 from the effective action as S e g(M) = — V^V^E, 0). Then the 
effective potential V^r is given by 



V cff (£,fl) = -ME cos 9 + In £ 

(l-4-2£ 2 sin 2 #) 1/2 



+ 



In 



2 />\ 1/2 



1+ (l -4-2E 2 sin 2 fl) 



(72) 



Now let us look into the vacuum structure of this effective potential by solving the saddle- 
point condition, which are given by 



cH4ff(E,fl) 
<9E 

86 



-Mcos6 + - 



8S sin 2 6 



2 fl \l/2 



E sin 9 



M 



s 1+ (l -4-2£ 2 sin 2 #) 
8£cos# 
1 + (l - 4- 2£ 2 sin 







2 /A 1 / 2 



. 



(73) 
(74) 
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Here we find two types of solutions for these equations depending on whether 9 is zero or 
nonzero: For a trivial solution 9 = 0, we have E = 1/M. For 9 ^ 0, the stationary conditions 
are written as 



ME - cos 9 = , (75) 
8E 2 

1 + (1 - 4- 2E 2 sin 2 #) 
Then, we find a solution for 9 ^ as 



1 " - ~ ■ — ■ 2M i/2 = • (76) 




(77) 
(78) 

8 — M 2 

Now we need to figure out which solution is realized as the vacuum of the theory by com- 
paring the potentials for the two solutions. We easily show for M 2 < 4, 

V eS (l/M,0)-V eS (Z,9) >0 . (79) 

while V eg (l/M, 0) -V eS (E, 9) < for M 2 > 4. Thus the vacuum of the strong-coupling QCD 
with the Hoelbling-type staggered- Wilson fermion is given by the following: For M 2 > 4 or 
equivalently M > —4^/2 + 2, M < —4\/2 — 2, there is only the chiral condensate as 

^(xx> = £cos£ =i, (80) 
i(xie xX ) = Esm0 =0. (81) 

For M 2 < 4 or equivalently — 4a/2 — 2 < M < — 4\/2 + 2, there is pion condensate which 
breaks parity symmetry spontaneously. 

1 M 

— (xx) = Scos^ = ^ , (82) 

iV XAA/ 8-M 2 



1 _ \2(A-M 2 ) 

— (xie x x) = Zsm9 = ±-± s . (83) 

iV xA/ 8-M 2 

The sign of the pion condensate Eq. (I53"j) reflects the Z 2 parity symmetry of the theory. The 
critical mass parameter M c = — 4y / 2 ± 2 and the range for the Aoki phase — 4a/2 — 2 < 
M < — 4v^2 + 2 is consistent with those of the hopping parameter expansion shown below 
Eq. fl30l . These results strongly suggest the existence of the parity-broken phase in the 
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FIG. 8: The pion condensate undergoes second-order phase transition. 

lattice QCD although it is just a strong-coupling limit. Figure M shows that the phase 
transition is second-order. 

We can also derive mass spectrum of mesons by expanding the effective action Eq. ( 1671) up 
to the quadratic terms of the meson-excitation field U x = Ai x — hA x . Since we are interested 
in the chiral limit which is taken from the parity-symmetric phase to the critical line, we 
here concentrate on the pion mass in the parity-symmetric phase. For the parity-symmetric 
phase (M 2 > 4), the quadratic part of the effective action is given by 

S cS (M) - S cS (M) =Y t S& ) (x,y)IIJL v 

^£(S n( - rtOT(p) - (84) 



where II (p) is the Fourier component of Il x , and 



v 1 



4 ^ 24 



(85) 



2S 2 

with p^± v = Pfj,±p v . Then we obtain the pion mass by solving V = at p = (im n a+7i, ir, n, n 



as 

2M 2 - 8 

cosh(m 7r a) = 1 H . (86) 

o 

By using the definition K = 1/2M with M = M + 2r and r = 2a/2, we find cosh(m 7r a) = 
1 + (1 — 16-fT 2 )/6A' 2 , which is consistent with the result of the hopping parameter expansion 
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Eq. ( 1351) : The pion mass becomes zero at the critical mass M 2 = 4, which indicates there 
occurs a second-order phase transition between parity-symmetric and broken phases in the 
strong-coupling limit. By defining quark mass as m q a = M — M c , we find PCAC relation 
near the critical mass as 

16 

(m n a) 2 = —m q a + G(a 2 ) . (87) 
3 

We can also study a case for non-zero spacial momenta by considering p = (iEa + ir,pia + 
7r,p20 + tt?P30 + in Eq. ( |85|) . By using the pion mass Eq. ( 1H71) and re- normalizing the 
Dirac operator as — |"D — > V, we show that Eq. fl85l) results in the Lorentz-covariant form 
up to 0(a) discretization errors, 

V = (E 2 — p 2 — ml)a 2 + G(a 3 ) , (88) 

with p 2 = p\ + p\ + p 2 . As we discussed in Sec. IHH we expected that we may find a 
sign of Lorentz symmetry breaking of the Hoelbling fermion in this study. However, we 
eventually cannot find a disease due to the symmetry breaking in the strong-coupling study. 
We consider that it is because Lorentz symmetry breaking appears mainly in the gluon 



sector as shown in [28| and it is difficult to find it in the meson sector in this limit. In future 
work, we may be able to find it by including higher order corrections of 1/g 2 and 1/N. 

We here discuss possibility of other condensation. For this purpose, we consider a general 
form of the meson field as 

M x = a + ie x 7T + £(-1)*^ + M-l^S + ^(-1)^ V > (89) 

where we define the vector, axial-vector and tensor meson fields as v^, a M and t M . We can 
easily show there is no other condensate by substituting this general form Eq. ( 1891) into the 
meson action Eq. fIBTj) . Thus we conclude that the vacuum we obtained is a true one. 

The results in this section suggest that the chiral limit can be taken in Hoelbling-fermion 
lattice QCD in a parallel manner to Wilson fermion. However we probably need to tune 
other parameters to restore Lorentz symmetry in lattice QCD with this type. Therefore, 
what we can state here is just that, if we succeed to restore Lorentz symmetry by parameter 
tuning, this fermion could be applied to lattice QCD as Wilson fermion. 
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B. Adams type 

We next investigate the case for the Adams fermion. We again consider the effective 
potential up to 0(Ai 3 ). The derivation is almost the same as the Hoelbling case in Sub- 
sec. IV Al The main difference between them is the number of the multi-links. The fermion 
of the Adams type includes the four-hopping terms while the Hoelbling one has the two- 
hopping terms. In the appendix |Bl we derive the effective potential for the Adams-type 
fermion. Here we only summarize the results. 

In this case, we again set r = 16v^3 to match the result to that of the hopping parameter 
expansion in Sec. IIVI We can derive the effective potential for scalar and pseudo-scalar fields 
by assuming condensation as M. x = 'Ee 16 * 6 . We note that the functional form of the effective 
potential is the same as Eq. f J72|) . By solving saddle-point equations, we find that the critical 
mass is given by M 2 C = 4 or equivalently M c = -16\/3 ± 2 with M = M + r and r = 16\/3. 
The vacuum in this case is given by following: For M 2 > 4 or M > —16-^/3 + 2, M < 
— 16-^/3 — 2, there is only the chiral condensate as 

±(x X ) = Zcos9 =-1, (90) 
^(xie x x) = Esin0 =0. (91) 

For M 2 < 4 or -16^ — 2 < M < -16^3 + 2, there emerge the pion condensate which 
breaks the parity symmetry spontaneously. 

1 M 

— (XX) = Scostf = ^ , (92) 

N x 1 8-M 2 



1 - - J2(A-M 2 ) 

-<fe> = £sin£ =±V__. (93) 

We note the critical mass and the parameter range of the Aoki phase are consistent with 
those of the hopping parameter expansion shown below Eq. (1401) . This result supports the 
existence of the parity-broken phase in the lattice QCD with the Adams fermion again. The 
behavior of pion condensate in this case is also given by Fig. [HI which shows the order of 
phase transition is second-order. This is consistent with the second-order scenario that we 
can take a chiral limit by a mass parameter tuning. 

We derive the pion mass from the quadratic parts of the effective potential in a parallel 
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way to the Hoelbling type. The pion mass for this case is given by 

v 2M 2 -8 

cosh(m 7r a) = H , (94) 

o 

for the parity-symmetric phase (M 2 > 4). By using the definition K = 1/2M with M = 
M + r and r = 16a/3, we find cosh(m 7r a) = 1 + (1 — 16-fT 2 )/10A' 2 which is consistent with 
the result of the hopping parameter expansion Eq. fj4"5]) : The pion mass becomes zero at the 
critical mass M 2 = 4, which indicates there occurs a second-order phase transition between 
parity-symmetric and broken phases in the strong-coupling limit. The PCAC relation holds 
near the critical mass also in this case. 

16 

(m^a) 2 = —m q a + 0(a 2 ) . (95) 



We can also show that the Lorentz-covariant dispersion relation recovers in the continuum 
limit in the Adams-type formalism as T> = (E 2 — p 2 — m 2 )a 2 + 0(a 3 ). It is reasonable that 
Lorentz-symmetric dispersion recovers since Adams fermion has sufficiently large discrete 
symmetry for Lorentz symmetry restoration. We can also argue the possibility of other 
condensations in the same way: We can show there is no other condensates by substituting 
a general form of the meson fields ( 189]) into the mesonic action for this case. 

All these results indicate that, in Adams- type staggered- Wilson, we can take a chiral limit 
by tuning a mass parameter toward the second-order critical line from the parity-symmetric 
phase. Since Adams fermion has sufficient discrete symmetry, it can be straightforwardly 
applied to lattice QCD in a parallel manner to Wilson fermions. 



VI. DISCUSSION ON TWO-FLAVOR CASE 



In this section, we discuss parity-flavor breaking for two-flavor staggered- Wilson fermions. 
We first consider two-flavor Hoelbling-type fermion action. In this case, except that the low 
discrete symmetry would require further parameter tuning, the situation is quite similar to 



that of the Wilson fermion 



12]. We here assume that mass parameters for two flavors Mf 



(f = 1,2) are equal, which means there is exact SU(2) flavor symmetry. The chiral and 
pion condensates are given by 

^ (XfXf) = £/ cos 9 f = 4- , (96) 

^ (XfiScXf) = s / sin d f = > ( 97 ) 
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Hoelbling type (two flavor) 



M c M 



Adams type 





(xiexX) = 








ml 







M c M 



FIG. 9: Conjectures of pion mass behaviors as a function of a mass parameter in two-flavor 
Hoelbling fermions and a single Adams fermion. In both cases PCAC relation holds in the parity 
symmetric phase. 

for Mf > 4 (parity-symmetric phase), while they are given by 



1 M f 

— (x f Xf) = Z f cos6f = is-, (98) 



1 - - \/2(4-M / 2 ) 

„<Xfie*Xf) = S/sin^ = ± V g _ Af2 / , (99) 

for Mj < 4 (Aoki phase). Here we have assumed that only the diagonal condensates in 
the flavor space (i.e. neutral condensates) can take finite values. Here we remind you of 
Mf = Mf + 2r with r = 2y2. We first look into the parity-symmetric phase. Although 
SU(2) chiral symmetry is explicitly broken due to the flavored-mass term, three-massless 
pions appear on the second-order phase boundary due to divergence of correlation length 
as shown in Fig. |9j In the parity-broken phase, things depend on whether or not 6\ and 82 
have the same sign in Eq. f )99|) . For Q\ =82, 



(xitxTiX) = , (i = 1,2,3) , 



(100) 
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where r, is the Pauli matrix and x stands a doublet x — (Xi? X2) T ■ For d\ = —9 2 , 

(xitxTix) = , 
(xitxT2X) = , 

(Wsx) 7^ . (101) 



rom Vafa-Witten's theorem 



30|, we expect that the latter vacuum (9i = —6 2 ) realizes 
16|. It is also possible to check this by studying next-leading-order calculation of 1/N 
or 1/g 2 expansions. If the latter scenario realizes, the flavored pion condensate Eq. 




breaks SU(2) flavor symmetry into its U(l) subgroup as well as parity symmetry 
Thus, in the parity-broken phase, we have two-massless pions as NG bosons associated with 
spontaneous breaking of the flavor symmetry. We summarize them in Fig. [9j This situation 



is qualitatively the same as the case of Wilson fermion 12J except possibility of further 
parameter tuning to recover Lorentz symmetry. 

The Adams-type staggered- Wilson fermion is more fascinating. It has two flavors for 
each branch in the first place. In this case there is no exact SU{2) flavor symmetry due 
to taste-mixing in original staggered fermions. It means that in the parity-broken phase 
there is no massless excitation since there is no continuous symmetry to be broken as Fig. |9j 
However the number of massless pions in the chiral limit (on the boundary) depends on 
residual discrete flavor symmetry in the pion sector. If the discrete flavor symmetry is not 
sufficient to have a degenerate pion triplet, we have only one massless pion in the chiral limit 
although these three are expected to be degenerate in the continuum limit. If the symmetry 
is large enough, we have three-massless pions on the phase boundary. The latter is a quite 
fascinating scenario because we can simulate two-flavor QCD with a single lattice fermion. 



It is possible to study it 
potential. Recently Ref. 



looking into transfer matrix symmetry or chiral Lagrangian 



28j | has reported that classification of pion operators from the 
transfer matrix symmetry indicates three-degenerate pions. Adams fermion can be a new 
standard of lattice fermion in the near future. 
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VII. SUMMARY AND DISCUSSION 

In this paper we investigate strong-coupling lattice QCD with staggered- Wilson fermions, 
with emphasis on the parity-broken phase (Aoki phase) structure. We consider hopping 
parameter expansion and effective potential analysis in the strong-coupling limit. We have 
shown that the parity-broken phase and the second-order phase boundary exist for both 
Adams-type and Hoelbling-type staggered- Wilson fermions, which is consistent with the 
second-order scenario for a chiral limit. 

In Sec. IIII| we discuss and classify discrete symmetries of two types of staggered- Wilson 
fermions. We show that they are invariant under charge conjugation and parity transfor- 
mation, the latter of which is defined as 4th-shift followed by spatial axis reversal. We also 
discuss smaller rotation symm etry in the Hoelbling fermion, which would require further 
parameter tuning as shown in 



In Sec. IIVt we analyze staggered- Wilson fermions by using hopping parameter expansion. 
From one-point functions of meson fields in the expansion, we find that pion condensate 
becomes nonzero in some range of the hopping parameter. From two-point functions, we 
show that square pion mass becomes zero on the boundary and becomes negative in the 
parameter region with nonzero pion condensate. These results suggest that there is a parity- 
broken phase and a second-order phase boundary. 

In Sec. |V] we study the effective potential for meson fields in the strong-coupling limit 
and large N limit to elaborate the phase structure in details. Here we develop a method to 
derive effective potential for lattice fermion actions with multiple-hopping terms. The gap 
equations from the effective potential exhibit nonzero pion condensate in the same parameter 
range as the hopping parameter expansion. From this analysis, we also show that pion 
becomes massless on the second-order phase boundary, and PCAC relation is reproduced 
around the boundary. If this property carries over into the weak-coupling regime, we can take 
a chiral limit by tuning a mass parameter in lattice QCD with staggered- Wilson fermions 
as with Wilson fermion. 

In Sec. IVIl we discuss the two-flavor cases. The situation in two-flavor Hoelbling-type 
fermions is similar to the original Wilson fermion except less rotational symmetry: Three- 
massless pions are expected to appear on the second-order critical lines, while two of them 
remain massless in the Aoki phase due to the flavored pion condensate. However we probably 
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need to care about Lorentz symmetry breaking in this case, thus we cannot straightforwardly 
apply it to two-flavor lattice QCD. The Adams-type staggered- Wilson fermion contains two 
flavors in each branch. Although the taste-mixing breaks flavor symmetry at finite lattice 
spacing, it does not necessarily mean non-degenerate three pions. Moreover SU(2) flavor 
symmetry should recover in the continuum limit at least, and three-massless pions emerge if 
we take a chiral and continuum limit. In this case, there is no rotational symmetry breaking, 
and the hypercubic symmetry will recover in the continuum limit. We can thus perform two- 
flavor QCD simulations with Adams-type staggered- Wilson fermion more efficiently than 
usual. 

All of these results shows new possibilities of lattice fermion formulations. In particular, 
the Adams fermion can be straightforwardly applied to 2-flavor lattice QCD since it does 
not require any other fine-tuning and automatically has two flavors. Taking account of less 
numerical expenses in the staggered fermion, there is possibilit y th at it would be numerically 
better than Wilson fermions, especially as an overlap kernel [H|. We finally note that the 
analysis here does not include contribution from some of higher-hopping terms or higher- 
meson fields. To confirm our results, we need to perform the same analysis with these higher 
contributions. In the future work, we can also study detailed mass spectrum of mesons and 
possibility of small other condensation in the Aoki phase. 
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Appendix A: spin and flavor separation 

From one-staggered field we define 16 species fields in the momentum space as <P{p)a = 
xip + t^a) i.— 71 /^ < Pp. < 7r /2) where tta (A = 1,2, ...,16) being 4-dim vectors whose 
components take or it. For convenience, we here consider a 16-multiplet field as <fi{p) = 
0(p) 2 , • • • ,<fi(p)w) T - As this 16-multiplet field has both the spinor (space-time) and 
the flavor (taste) indices, we can construct two sets of Clifford generators and S M , which 
operate on spinor and flavor spaces in the momentum field <fr(p). They satisfy the Clifford 
algebra as 

{r M ,r v } = 2^, (Ai) 

{H M ,H„} = 2(J Att/ , (A2) 
{r M ,5„} = 0. (A3) 

By using these definitions, the Dirac operator for the staggered fermion is given by D st 
iT^smpfj, for the 16 multiplet (j){p) |< 



Appendix B: Strong-coupling analysis for Adams-type 

In this chapter, we show the derivation of the effective potential for the Adams-type 
fermion in the strong-coupling limit. To derive the effective potential for the Adams type, 
we replace Eq. (|49|) by Eq. (1B1I) in the Adams type. 



exp 



U^+^o+p] exp [- (Tr{VEfl) - Tr(V^E))] . 

J \p^p^a J 

(Bl) 

Here we represent the partition function as 4 link integrals with U^x, U UtX +p, U PtX +p+i>, 
x+fi+o+p- V and E in Eq. ( 1B1I ) are the matrices which include components corresponding 
to 1-, 2-, 3-, and 4-link terms. The components of V and E consist of link variables and 
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TABLE III: Concrete forms of W ( jjJJ 9l ^ in Fig. H 
U ] M ] - -t/ f - . -fj 1 " * * * . 

p,x—p v,x—p.—v p,x—p,—v—p a,x—p—v—p—a 
^u,x— pt^ v i x ~ pVpiX— ji+oU(j,x— p+O+p 

+ u ux u ] , - -U ] . ~U ] . - . 

' fi^ v,x+p— v p,x+p—v—p a,x+p,—i/—p—a 

H I- + Up ]X U u ^~_~U p ^ x+ jx-c,U a ^ x j r jx-pj rp - 
I 77 1 77 ,r/t r/t 

T u p,x-p Uu < x ~P- u p,x-p+u-p u a,x-p+u-p-a 

+ + ■ ■ + U^ x U v>x jrpU px+ j 1+ ~_-U a)X +{ ii -i r p-p- 

I xft r/t 77 . .r/t 

p,x—p v,x — p—v p,x—p—u cr,x — p~u+p—a 

+ + H U^ !X U U:X+ fiUp !X +f l , + £,ul^ x+ j 1+ ~ + ~_ a 

+ [/ f -C/ 1 " - ~U ] - - -E^x-A-p-a 

p,x—p u,x—p—u p,x—p—u—p "l* P u P 

+ "I ^^^^^+^p,x+p+i>-p^i,x+p+u-p-a 

^ ' ^^^^x+^-u^P^+fi-^^x+fi-u+p-a 

~l I" Up ;X Ul x+ j 1 _ Up X+ j 1 _~_~U atX+ f l ~C'-p 

I" "I ^ ^,x-(P v ,^-P^ P,^-^+^ (T,x-jl+v+p-a 
I" + ^ p,,x-pPv,x-[l-v^ p,x-p\-uUcT,x-p,-i>+p 
~l I" ^\,x-p^^^-{^\,x~p^v~p^<y^-p^v~p 

fermion fields respectively. The concrete forms of the V and i? for this case are given by 

V$ = diag ■ • • , Vtf ) , (B2) 



diag {U^ x ) , 

diag (U ^^Uy^x+fiU p^x+ji+oUa^x+ji+o+f)) j 
diag {U^ X U V ^ X U P: x+i>Ua,x+o+p} , 
diag (u ^ x +vUl x+ jJJ p^+jiUa^x+fi+f^ ? 
diag (u ^x+pU u ,x+{i+pUl x+ii+0Ua,x+[i+i^ , 
diag (^U^ x +aU U) x+fj,+aU P: x+p+u+&Ulx + f l+0+ ^j 

diag ^^, x +/3+a^,s+A+p+CT^p, a; +/i+p+a^ r CT,x+/i+ 
diag i^^ x+ p + ^Ul x+ ^ +& U PtX+ ^ +Sr Ul x+jl+ A , 

diag (yl 

diag (ul x+i) U u , x+ pUl 



E% = diag (Ef, Ef, ■ ■ • , Eg) 



E x = diag (D ltli ) , 

^ = diag (D it p Upa ) , (z = 2, 3, • - ■ ,11) , 
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where we define the operator D as the fermion bilinears, 



~~ c^ln^XxXx+p i — 2 r lt 1 > x Xx+jiXx 5 



2, pupa 



D 



3, pupa 



^ A, pupa 
5, pupa 



ab 



ab 



) a ^ b ab b 

— ~S u X x +oXx+fi+p+a i — SuXx+fi+p+aXx+u 

b ab t 

= ~ s pX x +pXx+p+u+a ' (-^5,/Ji/) = S pXx+p+u+aXx+p 



S XxXx+p+u+p+a i (D2,pu) — SXx+p+u+p+aXx ) 



~ S P-Xx+pXx+O+p+a 



b,pupa 



pupa 



D 



\,li,v per 



D 



t 



9, pupa 



^10, pupa 



ab 



ab 



ab 



ab 



) a ^ b ab b 

~ s p+uX x +p+uX x +p+a 5 (-^10,/i!/) = s /i+i>X x + / 3+ ( 3-Xx+/iH-£' ' 



S <?X x +crX 



x+p+u+p ' 



(Ai, 



\ ab 



pu J 



s*X. 



x+p+u+p 



bXx+a 



~ s p+aXx+p+aXx+p+u 5 {D-j 



\ ab 



' p+aXx+p+uXx+p+a 1 



'Su+aXx+u+aXx+p+p i (D$,pv) — Su+aXx+p+pXx+O+a 5 
~ S u+pX x +u+pXx+p+a > (-D9,pu) — Sfr+pXx+p+aXx+v+p ' 



t \ Q '' b ab 

^11, pupa J = ~ S P+pXx+p+pXx+u+a > (Al1,|M/) = ^A+P^x+P+ctX 



(B17) 
(B18) 
(B19) 
(B20) 
(B21) 
(B22) 
(B23) 
(B24) 
(B25) 
(B26) 
(B27) 



Note that a = r (e m ) x /(4! • 16), s, = r (er) 5 ) x+fl /(4! • 16), and s, +u = r (er) 5 ) x+fi+i> /(4! ■ 16). 
Here V\ and E\ are 4x4 diagonal matrices while Vi and Ei (i — 2, 3, • ■ ■ , 11) are 24 x 24 
diagonal matrices. Here we note the situation of the cancellation between the diagrams 
crossing the different blocks is basically the same as the case for the Hoelbling type although 
there is difference between the 2-link and 4-link hoppings. We can derive W as a function 
of A by using the Schwinger-Dyson equation in a similar way to the Hoelbling type. A is 



A x = — 
16 



Y M x Mx + p + \Y M x+fl Mx 



+p+u 



p^u 



I 1 

- Y Mx+p+uMx+p+u+p + - Y Mx+p+o+pM 



x+p+0+p+cr 



p^u^p 

r 

4! • 16 



Y &M X M x+p+u+p+a + AMx+pMx+u+p+a + 2M x +p + uMx + p+a) ■ 



(B28) 
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Appendix C: Order parameter and zero eigenvalue of staggered- Wilson operator 



We investigate the relation between the order parameter (x^x^x) f° r spontaneous sym- 
metry breaking and zero eigenvalue of staggered- Wilson operator. In QCD, the chiral con- 
densate (ipip) which is the order parameter for spontaneous breaking of chiral symmetry 
relates to the zero eigenvalue of Dirac operator. It is called Banks-Casher relation 31] . In 
Wilson fermion, the pion condensate (ipi'j^T^ip) which is the order parameter for sponta- 
neous breaking of parity-flavor symmetry relates to the zero eigenvalue of Wilson operator. 
In staggered- Wilson fermion (two-flavor Hoelbling type), we derive the relation between 
(xie x T3x) and zero eigenvalue as Wilson fermion. Then we add the external field H for 
order parameter to the Hoelbling-type staggered- Wilson action Eq. ([12]) . 



S H (H) = xlDsw(M)+te x T 3 H] X 



(CI) 



The order parameter {x^xT^x) is represented as, 

lim (\ie x T-iY) = — lim lim — Tr [ ie x T-x 

-ff— >o 3A/ h->ov^ooV V D sw + ie x r 3 H 



— lim lim — tr 



it*. 



1 



— lim lim — tr 

V^-oo V 



Dsw 
1 



ie~H 



Dsw — i^xH 
1 



—% lim lim 

H^O V->oo V 



-i lim lim — 

V-s-oo V 



H SW + W 
1 



Hsw 
1 



iH 



dXp (A) 



X + iH X-iH 
1 1 



A) 



X + iH X-iH 



27rp(0) 



V 

= — lim lim lim 

e-S-0 H^fO V-»oo 



V 



(C2) 



where Hsw — txDsw is the Hermitian operator. Tr means the traces for flavor, color and 
space while tr means the traces for color and space. A and | A) are the eigenvalues and 
eigenstates and p(A) = J2y ^(^ ~ A') is the density of the state. By this analysis, we find 
the order parameter (x^x'TzX) f° r spontaneous breaking of parity-flavor symmetry relates to 
the zero eigenvalue of the staggered- Wilson operator Hsw- Also, we can derive this relation 
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for Adams fermion in the same way. 
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